Gaussian superpositions in scalar-tensor quantum 

cosmological models 



R. Colistete Jr.,* 

Laboratoire de Gravitation et Cosmologie Relativistes, Universite Pierre 
et Marie Curie, Tour 22, 4eme etage, Boite 142, 4 place Jussieu, 
p^. 75252 Paris Cedex 05, France 

J. C. FabW 

Departamento de Ffsica, Universidade Federal do Espfrito Santo, 
29060-900 - Vitoria, Espfrito Santo - Brazil 



m ; and N. Pinto-Neto* 

O ■ Centro Brasileiro de Pesquisas Ffsicas - Lafex, Rua Dr. Xavier 

in 
o 
o 
o 

7' 

u 

W). 

Abstract 



Sigaud 150, Urea 22290-180 - Rio de Janeiro, RJ - Brazil 

February 7, 2008 



A free scalar field minimally coupled to gravity model is quantized and the 
Wheeler-DeWitt equation in minisuperspace is solved analytically, exhibiting posi- 
tive and negative frequency modes. The analysis is performed for positive, negative 
and zero values of the curvature of the spatial section. Gaussian superpositions of 
the modes are constructed, and the quantum bohmian trajectories are determined 
in the framework of the Bohm-de Broglie interpretation of quantum cosmology. Os- 
cillating universes appear in all cases, but with a characteristic scale of the order 
of the Planck scale. Bouncing regular solutions emerge for the flat curvature case. 
They contract classically from infinity until a minimum size, where quantum effects 
become important acting as repulsive forces avoiding the singularity and creating an 
inflationary phase, expanding afterwards to an infinite size, approaching the classic- 
al expansion as long as the scale factor increases. These are non-singular solutions 
which are viable models to describe the early Universe. 
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1 Introduction 



The existence of an initial singularity is one of the major drawbacks of classical cosmology. 
In spite of the fact that the standard cosmological model, based in the classical general rel- 
ativity theory, has been successfully tested until the nucleosynthesis era (around t ~ Is), 
the extrapolation of this model to higher energies leads to a breakdown of the geometry 
in a finite cosmic time. This breakdown of the geometry may indicate that the classical 
theory must be replaced by a quantum theory of gravitation: quantum effects may avoid 
the presence of the singularity, leading to a complete regular cosmological model. 

The quantization of gravity is plagued with many conceptual and technical problems, 
and when it is applied to the whole universe new issues appear. In the Dirac quantization 
approach, a functional equation for the wave function of the Universe is obtained, the 
Wheeler-DeWitt equation 0, which is the basic equation of quantum cosmology. It is 
formulated in the so-called superspace, the space of all possible three-dimensional spatial 
geometries. It is very hard to find exact solutions of the full Wheeler-DeWitt equation, 
but solutions may be found in minisuperspaces where all but a finite number of degrees 
of freedom are frozen. 

Among the fundamental questions that come from the quantization of the universe as a 
whole, one of the most important concerns the interpretation of the wave function coming 
from the Wheeler-DeWitt equation. In order to extract predictions from the wave function 
of the Universe, the Bohm-de Broglie ontological interpretation of quantum mechanics 
0. |3|] has been proposed [|], |], 0> since it avoids many conceptual difficulties that follow 
from the application of the standard Copenhagen interpretation to an unique system 
that contains everything. In opposition to the latter one, the ontological interpretation 
does not need a classical domain outside the quantized system to generate the physical 
facts out of potentialities (the facts are there ab initio), and hence it can be applied to 
the universe as a whole[]With this interpretation in hands, one can ask if the quantum 
scenario predicted by the Wheeler-DeWitt equation is free of singularities and which type 
of classical universe emerges from the quantum phase. 

In a preceding work ||, we have applied this proposal to a free scalar-tensor model 
with minimal coupling in Friedmann-Robertson- Walker geometry, which can be obtained 
from a non-minimal scalar-tensor theory through a conformal transformation. Free scalar 
fields are good candidates to describe the material content of the early Universe because 
of their simplicity and because they represent stiff matter, the type of fluid advocated by 
Zel'dovich || to be relevant at early stages of cosmic evolution. Only positive curvature 
spatial sections have been studied. The bohmian trajectories in configuration space re- 
vealed an unexpected scenario: they behaved as the classical solutions for small values 
of the scale factor, but display quantum behaviour when the scale factor is big. As a 
consequence, the initial singularity is still present in this quantum model. 

The Wheeler-DeWitt solutions for this scalar-tensor model contain positive and neg- 
ative frequency modes, the first leading to an expanding universe, and the second to a 
contracting universe, near the singularity. Inspired by this observation, we constructed in 

1 Other alternative interpretations can be used in quantum cosmology like the many worlds interpre- 
tation of quantum mechanics W 
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H some particular superpositions mixing negative and positive models. In this way, we 
found non-singular quantum solutions which were, however, of planckian size and hence 
they could not be a model for our real Universe. 

The aim of the present work is to explore further the possibilities of the minisuper space 
model of Reference |J. First, we will not restrict ourselves to positive curvature spatial 
sections and second, we will explore more suitable superpositions of negative and positive 
modes, namely, the gaussian superposition. For the case the spatial section is flat, it is 
possible to solve analytically the expressions for the phase of the wave function, and to 
reduce the equations for the bohmian trajectories to a dynamical system. The critical 
points are calculated, and they are identified as center or nodes points. This leads to the 
existence of three kind of scenarios: periodic solutions representing oscillating universes; 
bouncing universes; models with a big-bang followed by a big-crunch. The bouncing 
universes contract classically from infinity until a minimum size, where quantum effects 
become important acting as repulsive forces avoiding the singularity, expanding afterwards 
to an infinite size, approaching the classical expansion as long as the scale factor increases. 
These are non-singular solutions which are viable models to describe the Universe we live 
in. For closed and open spatial sections, all calculations must be performed numerically, 
and the trajectories obtained in the configuration space reveal again the presence of 
oscillating universes besides those with a big-bang followed by a big-crunch. In all three 
cases, the oscillating universes have a characteristic scale of the order of the Planck length, 
except for very special gaussians in the case of zero saptial curvature. Hence, the most 
interesting scenarios emerge from the flat case, where we have succeeded to obtain a viable 
non-singular model. 

The article is organized as follows. In section 2, we describe the classical model and 
the corresponding Wheeler-DeWitt equation in the minisuperspace. Section 3 is devoted 
to the study of the gaussian superposition of the quantum solutions found before, and 
their corresponding analysis. In section 4 we present our conclusions. 



2 The classical and quantum minisuperspace models 

Let us take the lagrangian 



R - w(j>. p (j>> p 



For w = — 1 we have effective string theory without the Kalb-Ramond field. For w = —3/2 
we have a conformally coupled scalar field. Performing the conformal transformation 
Qiiv = ^9iiv we obtain the following lagrangian: 



-9 



R-[uj + 



(2) 



where the bars have been omitted. We will define C w = (u + §), which we will consider, 
from now on, to be strictly positive in order not to violate any of the energy conditions, 
at least classically. 
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We will consider the Robertson- Walker metric 



ds 2 = -N 2 dt 2 + 



aitf 
l + !r 2 



[dr 2 + r 2 (d# 2 + sin 2 (fl)V)] 



(3) 



where the spatial curvature e takes the values 0, 1,-1. Inserting this line element into the 
lagrangian and using the units where H = c = 1, we obtain the following action: 



S 



3V r Na 



4tt/2 



N 2 a 2 



+ C U 



6N 2 



(4) 



where V is the total volume divided by a 3 of the spacelike hypersurfaces, which are 
supposed to be closed, and l v is the Planck length. V depends on the value of e and 
on the topology of the hypersurfaces. For e = 0, V can have any value because the 
fundamental polyhedra of e = hypersurfaces can have arbitrary size (see Ref. [K|). In 



the case of e = 1 and topology S 3 , V = 2n 2 . Defining (5 2 = -ffi, 
again the bars, the hamiltonian reads 



2 

H = n( — f3 2 — + (3 
\ 2a 



2 Pi 
2a 3 



'2p 



and omitting 



(5) 



where 



Pa 



aa 



f3 2 N 
a 3 (p 



{3 2 N 



(6) 
(7) 



Usually, the scale factor has dimensions of length because we use angular coordinates in 
closed spaces. Hence we will define a dimensionless scale factor a = a/ (3. In that case the 
hamiltonian becomes, omitting the tilde: 



H 



N 
2P 



a 



ea 



As (3 appears as an overall multiplicative constant in the hamiltonian, we can set it equal 
to one without any loss of generality, keeping in mind that the scale factor which appears 
in the metric is (3a, not a. We can further simplify the hamiltonian by defining a = ln(a) 
obtaining 



H 



N 



where 



2 exp(3a) 



Pa 



2 i 2 
Pa +PS 



e 3a a 



eexp(4a) 



N 



,3a. 



N 



(9) 



(10) 
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The momentum p^ is a constant of motion which we will call k. The classical solutions 
are, in the gauge N — 1: 



1) e = 

= ±a + ci , (12) 
where C\ is an integration constant. In term of cosmic time they read: 

a = e a = 3kt 1/3 , (13) 
= ^ + c 2 . (14) 

The solutions contract or expand forever from a singularity, depending on the sign of k, 
without any inflationary epoch. 

2) e = 1 

k 



(15) 



cosh(20 — c\) 

where C\ is an integration constant, and from the conservation of p^ we get 

k = e 3a <P . (16) 

The cosmic time dependence is complicated and we will not write it here. These solutions 
describe universes expanding from a singularity till a maximum size and contracting again 
to a big crunch. Near the singularity, these solutions behave as in the flat case. There is 
no inflation. 

3) e = -1 

a = ea = v-^ — rr ' ^ 

| sinn (20 — c\) \ 

where C\ is an integration constant, and again, from the conservation of p^ we get 

k = e 3a (p . (18) 

As before, the cosmic time dependence is complicated and we will not write it here. These 
solutions describe universes contracting forever to or expanding forever from a singularity. 
Near the singularity, these solutions behave as in the flat case. There is no inflation^. 



Let us now quantize the model. The Wheeler-DeWitt equation is obtained through 
the Dirac quantization procedure where the wave function must be annihilated by the 

2 In the case e = — 1 there are classical solutions with C w < 0. Qualitatively, they represent universes 
contracting from an infinite to a minimum size and then expanding again to infinity 
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operator version of the constraint in Eq. With the choice of factor ordering which 
makes it covariant through field redefinitions, it reads 

-T^ + W + a * = ■ (19) 

Employing the separation of variables method, we obtain the general solution 

*(a,(f))= J F(k)A k (a)B k (<j>)dk , (20) 
where k is a separation constant, 

Bk{4>) = biexp(ik(J)) + b 2 exp(-ik(J)) , (21) 



and for e = 

for e = 1 

and for e = — 1 



A k (a) = ai exp(ika) + a 2 exp(— ika) , (22) 
A k (a) = a 1 I ik/2 (e 2a /2) + a 2 K ik/2 (e 2a /2) , (23) 
A k (a) = ai J lk/2 (e 2a /2) + a 2 N lk/2 (e 2a /2) . (24) 



The functions J, N, J, K are Bessel and modified Bessel functions of first and second kind. 

The Bohm-de Broglie interpretation of homogeneous minisuperspace models goes as 
follows: in general, the minisuperspace Wheeler-De Witt equation is 

H(p a (t),q a (m(q) = ■ (25) 
Writing \1> = Rexp(iS/h), and substituting it into (p5|) , we obtain the following equation: 

2 faM df a di, +U{q>l) + Q{qfl) = ' (26) 

where the quantum potential is 

«*>=-^*<5s • (27) 

The Bohm-de Broglie interpretation applied to quantum cosmology states that the 
trajectories q a {t) are real, independently of any observations. Eq. (p6|) is the Hamilton- 
Jacobi equation for them, which is the classical one amended with a quantum potential 
term (0), responsible for the quantum effects. This suggests to define 

P a = ^ , (28) 
dq a 

where the momenta are related to the velocities in the usual way: 

N dt K J 
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To obtain the quantum trajectories we have to solve the following system of first order 
differential equations, called the guidance relations: 



dS(q a ] 



dq a J N~ 

In the present case of the hamiltonian (||), the quantum potential 

,3a rQ2 R Q2 R , 



(30) 



becomes 



Q(a,< 



e 
2R 



and the guidance relations ( p0|) read 



as; 

da 
dS 



da 2 



e 3a a 
'~N~ 



(31) 



N 



(32) 



(33) 



Eqs. ( j30|) are invariant under time reparametrization. Hence, even at the quantum 
level, different choices of N(t) yield the same spacetime geometry for a given non-classical 
solution q a (t). There is no problem of time in the Bohm-de Broglie interpretation of 
minisuperspace quantum cosmology^. Let us then apply this interpretation to our min- 
isuperspace models and choose the gauge N — 1. 

3 Bohm interpretation of gaussian superpositions 

We will now make gaussian superpositions of these solutions and interpret the results 
using the Bohm-de Broglie interpretation of quantum mechanics. We will begin by the 
case e = 0, which is simpler, and it is the one to which the others reduce when a — > — oo. 



3.1 Hypersurfaces with e = 

This case can be solved analytically. The function F(k) is 

{k-df- 



F(k) = exp 
We can study two types of wave function: 



(7^ 



and 



#i(a,0) = J F(k)B k (<f))[A k (a) + A^(a)]dk 



* 2 (a,0) = / F(k)A k (a)[B k (<f>) + B_ k (<f>)]dk 



(34) 

(35) 
(36) 



3 This is not the case, however, for the full superspace (see Reference pi 
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both with a2 = &2 = 0. We will restrict ourselves to \&i because it yields the most 
interesting results. The results coming from \I/ 2 can be obtained from the first by changing 
a with 0. 

Performing the integration in k we obtain for vE^: 



= o\fn\ exp 



{a + 6) 2 a 2 - 



exp[id(a + (ft)} + exp 



[a 



)V 2 



exp[— id(a — 0)] 



(37) 

In order to obtain the bohmian trajectories, we have to calculate the phase S of the above 
wave function and substitute it into the guidance formula ([32 H33]) , working in the gauge 
N = 1. These equations constitute a planar system which can be easily studied: 





0<r 2 sin(2da) + 2<isinh(<7 2 a;0) 




exp (3a) < 


2[cos(2da) + cosh(a 2 a0)] J 



n = 7 — . (38) 



— acr 2 sin(2<ia) + 2<icos(2<ia) + 2<icosh(cr 2 a^ 



exp(3a)< 2[cos(2da) + cosh(<7 2 a<; 



(39) 



The line a = divides configuration space in two symmetric regions. The line = 
contains all singular points of this system, which are nodes and centers. The nodes 
appear when the denominator of the above equations, which is proportional to the norm 
of the wave function, is zero. No trajectory can pass through these points. They happen 
when = and cos(da) = 0, or a = (2n + l)ir/2d, n an integer, with separation n/d. 
The center points appear when the numerators are zero. They are given by = and 
a = 2<i[cotan(<ia)]/<T 2 . They are intercalated with the node points, and their separations 
cannot exceed n/d. As | a | — ^ oo these points tend to rni/d. As one can see from the above 
system, the classical solutions (a(t) oc i 1 / 3 ) are recovered when | a |— > oo or | |— ► oo, 
the other being different from zero. 

A field plot of this planar system is shown in Figure 1, for a — d — 1 . We can see plenty 
of different possibilities, depending on the initial conditions. Near the center points we can 
have oscillating universes without singularities and with amplitude of oscillation of order 
10. For negative values of a, the universe arises classically from a singularity but quantum 
effects become important forcing it to recollapse to another singularity, recovering classical 
behaviour near it. For positive values of a, the universe contracts classically but when 
and a are small enough, quantum effects become important creating an inflationary phase 
which avoids the singularity. The universe contracts to a minimum size and after reaching 
this point it expands forever, recovering the classical limit when a becomes sufficiently 
large. These are models which can represent the early Universe. We can see that for a 
negative we have classical limit for small scale factor while for a positive we have classical 
limit for big scale factor. 



4 As discussed above, these amplitudes can be very large as long as d becomes very small because the 
separation of the center points are of the order of 1 /d 
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-1.5 -1 -0.5 0.5 1 1.5 2 

a 

Figure 1: Field plot of the system of planar equations ( pI8| - P§D for a — d — 1, which 
uses the Bohm-de Broglie interpretation with the wave function ^i, Equation (|37]). Each 
arrow of the vector field is shaded according to its true length, black representing short 
vectors and white, long ones. The four shades of gray show the regions where the vector 
field is pointing to northeast, northwest, southeast or southwest. The black curves are 
the nullcline curves that separate these regions. The white points are the centers points 
whose neighbourhoods have oscillating trajectories. The trajectories are the white curves 
with direction arrows. 



For the wave function \I/2, the analysis goes in the same way but we have to inter- 
change a with (f). In this case we also have periodic solutions but the others are universes 
arising classically from a singularity, experiencing quantum effects in the middle of their 
expansion, and recovering their classical behaviour for large values of a. There are no 
further possibilities. 

We will now pass to the cases with curved spatial sections. One can immediately 
notice an important difference. The case e = has a symmetry a — > —a which is present 
not only in the Wheeler-DeWitt equation (|19"D but also in the solution ([37]). The cases 
e ^ do not possess this symmetry (the potential ee 4a in the Wheeler-DeWitt equation 
breaks it), and one should not expect to obtain the a > part of the field plots in these 
cases from the a < part through a reflection, as in the case e = (see Figure 1). 



3.2 Hypersurfaces with e = 1 

The Wheeler-DeWitt equation ( |T9"1) for e = 1, in the case we neglect the is analogous 

to a stationary Schroedinger equation with E = and V = e ia . Hence, one should make 
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superpositions involving only the parts of Ak(ot) which goes to zero as a goes to infinity, 
which are the Bessel functions K ik /2(e 2a ^ 2 ). Consequently, we will take the following 
superposition: 



^3 fa, 





f (k-d) 2 l 


j exp 


a 2 



K lk/2 (e 2a /2)B k (<f))dk . 



(40) 



The limit a — > — oo does not give the preceding results for e = because the Bessel 
function K reduces in this limit to 



K lk/2 (e 2a /2) « - 



exp[ik(a - ln(2)]r(l - — ) - exp[-ik(a - ln(2)]r(l + —) 



(41) 
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Figure 2: Field plot of the numerical solution of the system of planar equations (f44|-|4"5|) 
using the Bohm-de Broglie interpretation with the wave function \l/3, Equation (f40"D, for 
a = d = 1. Each arrow of the vector field is shaded according to its true length, black 
representing short vectors and white, long ones. The four shades of gray show the regions 
where the vector field is pointing to northeast, northwest, southeast or southwest. The 
black curves are the nullcline curves that separate these regions. The white point is the 
center point whose neighbourhood has oscillating trajectories. The trajectories are the 
white curves with direction arrows. 
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and the presence of the Gamma functions spoils their similarity. 
This case must be studied numerically, and the transformation 



a 



a — <p 
d> + a 



eases this task. The guidance relations (^2|-[33|) become 



a 



-2 exp 
-2 exp 



-3(a' + <P') 



-3(a' + d') 



dS_ 

dS 
da' 



(42) 
(43) 



(44) 
(45) 



and Figure 2 shows the field plot of this transformed planar system, using a — d — 
1. There are periodic solutions without singularities which happen when the bohmian 
trajectories cross the lines (a' = —3.73, <fi' < —3.73), or (<// = —3.73, a' < —3.73), or 
equivalently (a —\ (j) | —3.73, a < —3.73). These oscillating trajectories can reach very 
negative values of a but their maximum size cannot exceed a = 0, or a ~ l p i. Another 
behaviour is related to the trajectories shown in Figure 2 which are exclusively in the light 
gray region. They begin classically from a singularity, expand to a maximum value of a, 
and then return classically to a singularity. Concluding, we have two types of trajectories 
in this case: one which is periodic due to quantum effects, and the other which exhibit 
the pattern of classical behaviour: expansion from a singularity until a maximum size 
followed by a contraction to a big crunch. The periodic solutions have maximum size 
around a = 0, or a w l p i and they cannot represent the Universe we live in. 

3.3 Hypersurfaces with e = — 1 

In this case we will choose as AJa) the combination 



A k (a) 



r(l + y)J Jfc/2 (e 2 72) + r(i - l ±) J_ jfc/2 (e 2 72) 



(46) 



in order to get rid of the Gamma functions and obtain the preceding results for e = 
when a is very negative because the Bessel function J reduces in this limit to 



Jik/2(e 2a /2) 



exp[ik(a — ln(2)] 



r(i + f; 



(47) 



Taking this choice of A k (a), Eq. (f46|), into the gaussian superposition 
# 4 (a» = / exp - ^ -!- A k (a)B k (4>)dk , 



(4£ 



the numerical calculations with respect to the e = — 1 case show that the behaviour for 
very negative values of a is similar to the e = case, as one can see by comparing Figure 
3 with Figure 1. As a increases the regions with oscillating universes are squeezed and 
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Figure 3: Field plot of the numerical solution of the system of planar equations (p2|-|33|) 
using the Bohm-de Broglie interpretation with the wave function ^4, Equation (f48[), for 
a = d = 1. Each arrow of the vector field is shaded according to its true length, black 
representing short vectors and white, long ones. The four shades of gray show the regions 
where the vector field is pointing to northeast, northwest, southeast or southwest. The 
black curves are the nullcline curves that separate these regions. The white points are 
the centers points whose neighbourhoods have oscillating trajectories. The trajectories 
are the white curves with direction arrows. 



their separation decrease monotonically. Like the e = case, there are periodic solutions 
without singularities and with amplitude of oscillation of order 1. The other behaviour is 
described by trajectories that arise classically from a singularity, experiment a quantum 
halt at some maximum value of the scale factor, and then classically contracts to a big- 
crunch, contrary to the classical solutions of Eq. (|TT| ) which contract forever to or expand 
forever from a singularity. 



4 Conclusion 

The quantization of a scalar-tensor model in the minisuperspace leads to a separable 
partial differential equation, admitting analytical solutions, with positive and negative 
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frequencies. In this work, we have studied gaussian superpositions of these different 
modes and the corresponding bohmian trajectories. Such analysis was performed for zero, 
positive and negative curvature spatial sections, which are considered to be compact. The 
bohmian trajectories in configuration space were calculated numerically, excepted for the 
flat case, where it is possible to reduce the equations for the bohmian trajectories to a 
two dimensional dynamical system. 

The comparison of the trajectories in the configuration space of the variables a and 
<ft, which are the dynamical degrees of freedom of the minisuperspace, with the classical 
ones, allows one to identify the classical and quantum phases for the scalar-tensor cos- 
mological models. For all three different values of the curvature of the spatial sections, 
the configuration space of the quantum solutions displays oscillating universes. However, 
these oscillating universes remain at the Planck scale and they cannot be considered as 
candidates for the description of the early Universe (they are more like baby universes), 
except for the unnatural choice \d\ << 1 in the e = case. There are also trajectories 
which correspond to universes which begin and end in singular states. Only for the flat 
case it is possible to have bouncing models. 

In the bouncing models of the flat spatial section case, the scale factor has an infinite 
initial and final values, near which it behaves classically. As it approaches the singularity, 
the repulsive quantum effects lead to the bounce, avoiding the singularity. Such a scenario 
can be a candidate for the description of our early Universe, since it is free from the initial 
singularity and behaves classically in the asymptotic limit of large values for the scale 
factor. It is worth to note that this classical asymptotic limit corresponds to the stiff 
matter Friedmann universe which, according to Zel'dovich ||, is the most promising one 
to describe the very early Universe. 

The free scalar field model considered in this paper, on the other hand, can be con- 
nected to a non- minimal coupled scalar field, with a coupling parameter u, like in the 
Brans-Dicke theory, by a conformal transformation. In [12| a quantum analysis of these 
non-minimal models was performed, and it was shown that non-singular scenarios can be 
obtained when the parameter u is negative. In fact, all quantum analysis performed here 
can be connected with a similar analysis in the non-minimal case through a conformal 
transformation, namely a^MC = a + and 4>nmc — 0- One can verify that when the 
minimal model displays singularities, it is possible to have non-singular solutions in the 
corresponding non-minimal case; but the non-singular solutions in the minimal case must 
also be non-singular in the corresponding non-minimal models. 

An important generalization of the model studied here would be to consider self- 
interacting scalar fields. It was shown in []nj that to each perfect fluid barotropic equation 
of state, it is possible to construct a self-interacting scalar field model leading to the same 
classical description in minisuperspace. We may argue if this correspondence remains 
at the quantum level. Note that in Ref. we have obtained bouncing universes for 
radiation fields with e = and e = — 1 with the same qualitative behaviour as the 
bouncing universes found here. They are also viable models for the early Universe. One 
should investigate if a scalar field model with a potential corresponding to the radiation 
fluid would give similar results. 
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